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Multivalley spin relaxation in n-type bulk GaAs in the presence of high electric field is investigated 
from the microscopic kinetic spin Bloch equation approach with the V and L valleys included. We 
show that the spin relaxation time decreases monotonically with the electric field, which differs from 
the two-dimensional case and is recognized due to to the cubic form of the Dresselhauss spin-orbit 
coupling of the V valley in bulk. In addition to the direct modulation of the spin relaxation time, the 
electric field also strongly influences the density and temperature dependences of the spin relaxation. 
In contrast to the monotonic decrease with increasing lattice temperature in the field-free condition, 
the spin relaxation time is shown to decrease more slowly under the influence of the electric field 
and even to increase monotonically in the case with small electron density and high electric field. 
We even predict a peak in weakly doped samples under moderate electric field due to the anomalous 
lattice-temperature dependence of the hot-electron temperature. As for the L valleys, we show that 
instead of playing the role of a "drain" of the total spin polarization as in quantum well systems, in 
bulk they serve as a "momentum damping area", which prevents electrons from drifting to higher 
momentum states. This tends to suppress the inhomogeneous broadening and hence leads to an 
increase of the spin relaxation time. 

PACS numbers: 72.25.Rb, 72.20.Ht, 71.10.-w, 71.70.Ej 



I. INTRODUCTION 

Semiconductor spintronics, which aims at incorporat- 
ing the spin degree of freedom in electronics, has at- 
tracted much attention in the past decades.— ~— One of 
the main challenges in realizing spintronic devices lies 
in the control of spin lifetime. Many investigations 
have been devoted in bulk III-V semiconductors 
Various factors on spin relaxation, such as carrier and 
impurity densities^— temperature ^r— mobility^ ini- 
tial spin polarization^ and electric 20 ' 28 and magnetic 
fields r°- ~ 12 d 7 i 24 etc, have been extensively studied. Re- 
markably, spin relaxation time (SRT) as long as 130 ns 
is observed in n-type GaA o 10 ' 11 and the SRT varying 
by more than three orders of magnitude with tempera- 
ture or density is reported^ Very recently, Jiang and 
Wu performed a systematic investigation on the elec- 
tron spin relaxation in bulk III-V semiconductors in 
metallic regime from the fully microscopic kinetic spin 
Bloch equation (KSBE) approach^ The D'yakonov- 
Perel' (DP) mechanism^! is demonstrated to be dom- 
inant in n-type III-V semiconductors, even in narrow 
band ones, and many features in contrast to the previ- 
ous understandings in the literature are reported^ Some 
of their predictions have been soon realized in the sub- 
sequent experimental investigations ^ 2 - ~ 16 ' 24 " — Moreover, 
the very recent works 29 '- on the spin relaxation near 
the metal-to-insulator transition have extended the in- 
vestigation in that regime back in 2002. 11 

It is noted that works discussed above mainly concern 
electron spin relaxation in the absence of or with rel- 
atively low electric fields. In Ref. [2(], an electric field 
applied up to 2 kV/cm is also considered and the SRT 



is found to be effectively manipulated, especially for the 
low temperature case, where the SRT is suppressed down 
to 1/10 of its original value with a small electric field up 
to 0.05 kV/cm. Clear hot-electron effect is shown and 
the different field dependence of the SRT compared to 
that in quantum wells is addressed. Nevertheless, since 
most current electronic devices are performed in the high- 
field conditions, new features of the SRT are expected 
when electrons are drifted to high valleys 22- The differ- 
ent spin-orbit couplings (SOCs), momentum relaxations 
and effective masses in different valleys should have pro- 
nounced effects on spin dynamics. In fact, investigations 
on how the spin relaxation is affected by high electric 
field have been carried out in quantum well systems. The 
hot-electron effect^ and the multi-subband effect^ on 
spin dephasing in n-type GaAs quantum wells have been 
studied. Different field dependences of SRT are observed 
for different temperatures, well widths and initial spin 
polarizations. Moreover, the multivalley spin relaxation 
under high in-plane electric fields has been investigated 
in n-type GaAs quantum wells by taking into account 
the r and L valleys/- 3 - It is predicted that although the 
SOC in the L valleys is much larger than that of the T 
valley,— the spin polarization in the L valleys shares the 
same damping rate as that in the T valley. The L valleys 
are pointed out to play the role of a "drain" of the total 
spin polarization due to the large SOC therein. Non- 
monotonic dependence of the SRT on electric field has 
been reported, while the spin Gunn effect^ the sponta- 
neous spin-polarization generation in the high-electric- 
field charge Gunn region, is pointed out to hardly re- 
alize in GaAs quantum wells. Despite these works in 
two-dimensional structures, a detailed fully microscopic 
study in bulk system under high electric field is still ab- 
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sent. How the temperature and electron density depen- 
dences of the SRT are affected by the high electric field 
is still unclear and the feasibility of the spin Gunn ef- 
fect predicted in bulk 3 - 7 - needs to be checked. Because of 
the cubic form of the Dresselhauss SOC in the T valley 
as well as the absence of quantum confinement in bulk, 
different behaviors in spin relaxation are expected. 

In this work, we apply the KSBEs to investigate the 
electron spin relaxation in n-type bulk GaAs in the pres- 
ence of high electric fields with the T and L valleys in- 
cluded. The electric field dependence of the SRT is calcu- 
lated and found to be monotonic, in contrast with the re- 
sults in quantum wells^ We attribute this to the strong 
enhancement of inhomogeneous broadening 7 - ^ 8 ' 39 result- 
ing from the cubic form of the T- valley Dresselhaus SOC 
in bulk. The electric field is also shown to effectively 
change the density and the temperature dependences of 
the SRT. Remarkably, a peak is predicted in the tem- 
perature dependence of the SRT for relative low electron 
density under a moderate electric field. The underlying 
physics is analyzed. As for the L valleys, we show that 
in contrast to the two-dimensional system where the L 
valleys play the role of a "drain" of the total spin po- 
larization, unexpectedly in bulk their effect on the total 
spin dynamics is weak and the total spin relaxation is 
mainly determined by the T valley. It is shown that even 
with much larger SOC in the L valleys, the inclusion of 
the L valleys results in a longer rather than a shorter 
SRT. We find that the L valleys serve as a "momen- 
tum damping area" where electrons are hardly to drift to 
higher momentum states due their large effective mass. 
This tends to suppress the inhomogeneous broadening 
and hence leads to the increase of the SRT. 

This paper is organized as follows: In Sec. II, we set 
up our model and construct the KSBEs. In Sec. Ill, we 
lay out our main numerical results. The effects of electric 
fields, together with the carrier density, temperature and 
intervalley scattering (note that in this paper, with the 
term "intervalley scattering" , we always indicate the T-L 
intervalley scattering unless otherwise specified) on the 
SRT are discussed. We summarize in Sec. IV. 



the L valleys, since even for the highest electric field 
(E = 8 kV/cm) employed in this work, the fractions of 
electrons in these valleys are negligible^ 

The KSBEs derived from the noncquilibrium Green 
function method read o 7 ' 35 ' 38 ' 39 
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in which pAk A is the density matrix of electrons with 
momentum k> in A valley. Note that h\ is defined in 
reference to the valley center in each valley. The diag- 
onal term pxk x ,aa = f\k x ,cr (c = ±1/2) describes the 
distribution of each spin band and the off-diagonal term 
P\\, i i = pX, i i is the correlation between the two 
spin bands. The coherent term is given by 



dtpxk x | 



^A(k> 



££ F (k A ) 



PAk A 



(2) 



where [ , ] is the commutator and SIa^a) represents the 
Dresselhaus SOC in A valley*^ 3 - By setting the [100] and 
[010] directions as x- and y-axises, respectively, f2r(kr) 
takes the form 43 ' 44 
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By further denoting = (1, 1, ±l)/y3 and n 2 /4 = 

-(I,±l,l)/\/3, which are the unit vectors of the longi- 
tudinal principle axis of the Li valleys, we have for the L 
valleys^ 3 -^ 
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~ Ek A ^AA : k A -k' A PAk' A in Eq. {2]) is the 
Coulomb Hatree-Fock term with V\A,k A -k A representing 
the intravalley Coulomb scattering matrix element. The 
drift term takes the form <9 t pAk A | drift = eE • V kA PAk A 
(e > 0). <9t/9Ak A | scat stands for the scattering terms, 
which include intra- and inter-valley parts with the ex- 
plicit expressions given in Appendix A. 



II. MODEL AND KSBES 

We start our investigation in n-type bulk GaAs where 
the four L valleys locate at L points [K°. = ^-(1, ±1, ±1) 
with ao denoting the lattice constant and i = 1-4] and 
lie energetically Etl — 0.296 eV above the T valley It 
is noted that the four L valleys can be arbitrarily cho- 
sen from the eight L points limited by the condition that 
there is no center inversion symmetry between any two 
of them. In the spherically symmetric approximation, 
the electron effective masses of the T and L valleys are 
represented as mr = 0.067mo and rriL — O^too ) 41 ' 42 re- 
spectively, with too representing the free electron mass. 
We do not consider valleys of higher energy, e.g., the 
next-nearest valleys X which are Elx = 0.166 eV above 



III. NUMERICAL RESULTS 

In this section, we present our results obtained by nu- 
merically solving the KSBEs following the scheme laid 
out in Refs. [20| and l35l . All parameters used in our com- 
putation are listed in Table I. 

The initial spin polarized state of the system is pre- 
pared in the similar way as in Ref. |3o1 . starting from the 
unpolarized equilibrium state. That state is first driven 
to the drifted steady state under the electric field. The 
main difference lies that after driving the system to the 
steady state, we turn on a circularly polarized laser pulse 
at t\ — 6 ps to excite spin polarized electrons with a de- 
gree of injected spin polarization Pinjoct = 50% into the 
r valley^ Sf rkr ,a=a a exp[-(el r - e pump ) 2 /2S 2 e ] [1 - 

/rk r ,<r(ii)]- Here a a = n pumPt(T /{J2k r e xp[-( £ k r ~ 
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2 /2<5 2 ][l — /rk r ,o-(*i)]}- £pum P is the energy of pulse 8 kV/cm where the X valleys can be neglected^ 



pump,2 

= 4 meV 



center in reference to the band minimum and 5 e =H/6 T 
with S T denoting the pulse width. iip Ump .„ is the density 
of electrons with spin a after excitation with n c 
3n pump _i. In this work, we employ £ pump 
for the case of resonant excitation, 5 T = 0.01 ps and 
n pump = n P ump,i + "pump -j = 0-02 x n e . Here n e 
stands for the unpolarized electron density before pump- 
ing which is equal to the doping density, so the total 
spin polarization after the pump pulse is Pq ~ 1%. It is 
noted that due to the strong electron-electron Coulomb 
scattering, the drifted steady-state distribution is estab- 
lished within 0.1 ps after the pump pulse. Therefore the 
pulse characters, i.e., £ pump and 5 T , have little influence 
on the SRT. It is further noted that due to the small ini- 
tial spin polarization, the exact value of the polarization 
has marginal effect on the SRT. 20 

In Fig.[TJ we plot the typical time evolution of the spin 
polarization and also the electron population in A valley 
in the condition with electric field E = 6 kV/cm, electron 
density n e = 10 16 cm -3 and temperature T = 300 K. 
Clear transfer of electrons from the T valley to the L val- 
leys is observed. In order to quantitatively understand 
the influence of the L valleys, we calculate the steady- 
state drift velocities of A-valley electrons, the mobilities, 
the electron populations and the hot-electron tempera- 
tures in the T and L valleys as function of electric field. 
The explicit results are plotted in Appendix B with un- 
derlying physics analyzed. From Fig.[5Ja) of Appendix B, 
one finds that our model is validated against the experi- 
mental data^ 
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FIG. 1: (Color online) Typical evolution of the electron spin 
polarization and population in A valley in the condition with 
n e = 10 16 cm" 3 , E = 6 kV/cm and T = 300 K. 

In the following, we focus on how the SRT is influenced 
by the electric field, carrier density and lattice temper- 
ature. The effect of the L valleys is also studied. The 
SRT t is obtained by fitting temporal evolution of spin 
polarization P\ with an exponential decay (see the inset 
of Fig. [TJ . Throughout this work, the electric field is ap- 
plied in the x-direction and its strength is limited up to 



TABLE I: Material parameters used in the calculation (from 
Refs. [3^ and [13 unless otherwise specified). 



(eV) 
(eV) 



Ef (eV) 
mf/mo 
m* L /m 

Kg 

r^OO 

7 d (eV-A 3 ) 
/3d (eV-A) 



1.519 

1.815 

1.981 

0.067 

0.23 

12.9 

10.8 

23.9 a 

0.26 



firr (meV) 
n LiLl (meV) 
QrL (meV) 
fi LiLj (meV) 
£>L iLl (10 9 eV/cm) 
D F l (10 9 eV/cm) 
D LiLj (10 9 eV/cm) 
d (10 3 kg/m 3 ) 



35.4 

34.3 

20.8 

29.0 

0.3 

1.1 

1.0 

5.36 



Ref. |36J. 



A. Electric field dependence and effect of L valleys 

We first study the electric field dependence of the SRT 
at T = 300 K. In Fig. El we plot the SRTs against 
electric field for electron densities n e = 10 16 , 10 17 and 
10 18 cm -3 , corresponding to the nondegenerate, interme- 
diate and degenerate regimes in the field-free condition, 
respectively. It is seen that the SRT is effectively modu- 
lated by the electric field. For all three cases, the SRTs 
decrease monotonically with the electric field and reach 
down to the ones with one order of magnitude smaller 
than the corresponding field- free values at E = 8 kV/cm. 
This is different from the previous work in n-type GaAs 
quantum wells where a nonmonotonic electric field de- 
pendence is observed^ 3 -"— but consistent with the very 
recent work in bulk GaAs by Jiang and Wu where the 
electric field is applied up to 2 kV/cm^ In the regime 
with E < 2 kV/cm, according to Fig.[Btc) in Appendix B, 
the L valleys are still irrelevant. Therefore the underly- 
ing physics of the monotonic decrease of the SRT is the 
same as that analyzed in Ref. Due to the cubic form 
of the Dresselhauss SOC in bulk, the enhancement of 
inhomogeneous broadening from the drift effect and the 
hot-electron effect [see Fig.[B]Jd) in Appendix B where the 
two-hot-electron-temperature behavior of T-valley elec- 
trons is discussed] is more pronounced than that in quan- 
tum wells with small well width where the SOC is in the 
linear formi^ 3 -"— It hence overtakes the enhancement of 
momentum scattering and leads to the monotonic field 
dependence of the SRT. However, when the electric field 
is increased over 2 kV/cm, a visible amount of electrons 
are drifted into the L valleys [see Fig.lBJc) in Appendix B] 
and hence the L valleys are expected to play a role in the 
total spin relaxation. 

An important feature of the multivalley spin relaxation 
can be noticed in the inset of Fig. [T] that the evolutions 
of spin polarizations in the T and L valleys share the 
same damping rate. This can be understood in respect 
to the rapid exchange of electrons between the F and 
L valleys resulting from the strong intervalley scatter- 
ing. As shown in Fig. [21a,), by removing the intervalley 
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FIG. 2: (Color online) SRT r vs. electric field E at temper- 
ature T = 300 K for three electron densities: n e = 10 16 , 10 17 
and 10 18 cm" 3 . 



scattering Hyl, one finds that the spin polarizations in 
separate valleys now evolve independently and the "in- 
trinsic" SRT in the L valleys (curves with o) becomes 
one order of magnitude smaller than that of the F val- 
ley (curve with □). This demonstrates the crucial role 
of the intervalley scattering in obtaining the identical 
damping of spin polarizations in the T and L valleys. 
To be more specific, we also calculate the SRT in each 
valley by removing only the intervalley electron-phonon 
or electron-electron Coulomb scattering. It is seen from 
Fig. EJa) that the SRTs in the T and L valleys with- 
out the intervalley electron-electron Coulomb scattering 
coincide, whereas those without the intervalley electron- 
phonon scattering are far away from each other. This 
demonstrates that the identical damping of spin polar- 
izations in the T and L valleys comes from the intervalley 
electron-phonon scattering. It is noted that the fast spin 
relaxation in the L valleys hints that they may serve as 
a "drain" of the total spin polarization, just as the case 
in quantum wells ^ 

To further elucidate the role of the L valleys in the bulk 
system, we calculate the SRTs by artificially removing 
the L valleys [labeled as Case (ii)] and by setting the SOC 
in the L valleys to zero [labeled as Case (iii)] in the com- 
putation. The calculated SRTs are compared to those 
obtained in the genuine condition [labeled as Case (i)]. 
In Fig. |31[b) , we plot the SRTs as function of electric field 
for Case (i)-(iii) with n e = 10 18 cm" 3 and T = 300 K. 
One observes that for E < 2 kV/cm, the effect of the L 
valleys is marginal due to little electrons drifted into the 
L valleys under weak electric field. When the electric field 
further increases, the L valleys start to play an important 
role and obvious distinctions are seen for the three cases. 
We first compare the genuine SRT with the one without 
the L valleys. As shown in Fig. [21b), by removing the 
L valleys in the calculation, surprisingly the SRT (curve 
with •) becomes smaller instead of larger than its gen- 
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FIG. 3: (Color online) For (a), (b), (d) and (e): in bulk 
with electron density n e = 10 18 cm -3 : (a) SRTs r of the V 
and L valleys as function of electric field E by only removing 
the intervalley scattering Hlf (curves with □ and o), the 
intervalley electron-phonon scattering i?2r (curves with A 
and •) or the intervalley electron-electron Coulomb scattering 
HTv (curves with V and x ) . (b) SRT r vs. electric field E in 
the genuine condition [Case (i), curve with □], by setting the 
SOC in the L valleys to zero [Case (iii), curve with A] and 
without the L valleys [Case (ii), curve with •]. (d) and (e) 
Mobility /i and hot-electron temperature Tr vs. electric field 
with [Case (i), curve with □] and without [Case (ii), curve 
with •] the L valleys. For (c) in quantum well (QW) with 
well width a = 6 nm and n e — 6 x 10 11 cm -2 , SRT r vs. 
electric field E in the genuine condition [Case (i), curve with 
□], and by setting the SOC in the L valleys to zero [Case (iii), 
curve with Al. T = 300 K in all the situations. 



uine value (curve with □). For E = 7 kV/cm, the SRT 
without the L valleys reaches down to half of the genuine 
value. In order to understand this behavior, we calculate 
the mobility and the hot-electron temperature Tr , which 
serve as scales of the drift and hot-electron effect, respec- 
tively. From Fig. E^d) and (e) , pronounced distinctions 
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are seen for E > 4 kV/cm. In the case without the L 
valleys, the mobility increases rather than decreases with 
the electric field and reaches more than twice of its gen- 
uine value at E = 7 kV/cm. Meanwhile, the hot-electron 
temperature Xr reaches almost twice of its genuine value. 
Consequently, both the drift and hot-electron effects are 
markedly enhanced without the L valleys, which in turn 
leads to a drastic increase of the inhomogeneous broaden- 
ing due to the cubic form the T-valley Dresselhaus SOC 
in bulk. It overtakes the effect of the absence of the L 
valleys and leads to the smaller SRT compared to the 
case with the L valleys. This tells us that the L valleys 
serve as a "momentum damping area" where electrons 
from the T valley are blocked from reaching higher mo- 
mentum states, thanks to the large effective mass in the 
L valleys. 

We then focus on the curves corresponding to Case 
(i) (curve with □) and (iii) (curve with A). The SRT 
calculated by setting the SOC in the L valleys to zero 
becomes a little larger and reaches less than 3/2 of the 
corresponding genuine value at E = 7 kV/cm, indicat- 
ing that the spin relaxation in the L valleys only slightly 
modulates the total spin relaxation. This is very differ- 
ent from the quantum well system. As a comparison, 
in Fig. [31(c) we plot the SRTs against electric field in a 
quantum well in Case (i) and (iii) i 35 i 46 The well width is 
chosen as 6 nm and the two-dimensional electron density 
is n e — 6 x 10 11 cm -2 . It is seen that the behavior of 
field dependence of the SRT is totally changed in the ab- 
sence of the SOC in the L valleys. Instead of decreasing 
rapidly after a small increase, the SRT with = in- 
creases monotonically and reaches more than ten times 
of the corresponding genuine value at E = 7 kV/cm. 

The underlying physics is understood as follows: in 
bulk, due to the cubic form of the Dresselhauss SOC 
in the T valley, the field-induced drift effect and the hot- 
electron effect lead to a stronger enhancement of the inho- 
mogeneous broadening compared to the increase of mo- 
mentum scattering. Moreover, a simple calculation shows 
that although in average, the inhomogeneous broadening 
in the T valley is much smaller than that in the L valleys, 
in the energy range of the T valley overlapping with the 
L valleys (which is roughly where the exchange of elec- 
trons between the T and L valleys happens) , the effective 
magnetic field from the Dresselhauss SOC is already com- 
parable with that in the L valleys. However, in quantum 
wells with a small well width, the enhancement of scat- 
tering is more effective than that of the inhomogeneous 
broadening thanks to the linear form of the Dresselhaus 
SOC. 55 Besides, the effective magnetic field in the energy 
range of the T valley overlapping with the L valleys is still 
much smaller than that in the L valleys. In addition, 
it is further noted that electrons are more liable to be 
drifted into the L valleys in quantum wells. In Fig. [31(c), 
rii/rie reaches 45% at E = 7 kV/cm compared to 25% in 
bulk, which further enhances the effect of the L valleys 
in quantum wells. Nevertheless, in Fig. EIc), the electric 
field corresponding to ni,/n e = 0.25 is E — 5 kV/cm, at 



which the SRT without the SOC in the L valleys is still 
about five times larger than the corresponding genuine 
value. All these lead to the pronounced difference be- 
tween two- and three-dimensional systems. It also tells 
us that the L valleys no longer serve as "spin drain" in 
bulk, differing from the case of quantum wells^ 5 - 



B. Density dependence 

Another interesting feature seen in Fig. [2] is that the 
t-E curve corresponding to the electron density n e — 
10 18 cm~ 3 intersects with the other two, indicating that 
the density dependence of the SRT changes with the vari- 
ation of electric field. In order to elucidate this behavior, 
we calculate the SRT by varying the electron density n e 
from 2 x 10 15 to I0 18 cm~ 3 with electric field E = 0, 4 
and 6 kV/cm. Note that in bulk the impurity density rii 
is always taken as n e in this paper. From Fig. 31 one ob- 
serves that for E = kV/cm, there is a peak at around 
n e — 2.0 x I0 17 cm~ 3 . This peak in the density depen- 
dence of the SRT has been theoretically predicted^ and 
experimentally confirmecU 3 - -15 ^ 4 very recently, and is at- 
tributed to the crossover from the nondegenerate to de- 
generate limit when the electron-impurity and electron- 
electron Coulomb scatterings are the dominant scattering 
processes^ The crossover between degenerate and non- 
degenerate limits can be estimated by the Fermi temper- 
ature Tp [= (37r 2 n e ) 2//3 / (2m)], with the peak determined 
by Tp ~ iT-TJ^r^ML-^. Here the Fermi tempera- 
ture corresponding to the electron density at the peak 
is Tp w 207 K, comparable to the lattice temperature 
T = 300 K (Tp /T ~ 2/3), in line with the previous 
works. 

By applying an electric field E = 4 kV/cm, one 
observes from Fig. [4] that the peak moves to around 
n e = 5 x I0 17 cm~ 3 , and when the electric field is further 
increased to E — 6 kV/cm, the peak is shifted to even 
higher electron density (at around n e = 9 x I0 17 cm~ 3 
hence is not very obvious in the figure). We point out 
that this is due to the increase of hot-electron temper- 
ature induced by the strengthened electric field. As 
pointed out by Shenjii because of the laser-induced hot- 
electron effect, the peak of density dependence of the SRT 
appears at where the hot-electron temperature T e , rather 
than the lattice temperature T, approximately equals the 
Fermi temperature Tp. The underlying physics is sim- 
ilar here with the electric-field-induced hot-electron ef- 
fect. However, as shown in Fig. [5] in Appendix B, the 
L valleys start to play a role with high electric field and 
moreover, the two-hot-electron-temperature behavior of 
the T-valley electrons is seen. This makes the situation 
more complex. Nevertheless, by comparing the curves of 
"intrinsic" SRTs of the T and L valleys in Fig. [31(a) with 
the corresponding genuine one in Fig. [2Kb), one notices 
that the curve of the "intrinsic" SRT in the T valley re- 
sembles the genuine one while that of the "intrinsic" SRT 
in the L valleys is far away it. This indicates that the 
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the multivalley spin relaxation in the presence of high 
electric field is mainly determined by the T valley. Fur- 
thermore, since most electrons stay in the T valley for 
electric field up to 6 kV/cm [see Fig.^c) in Appendix B] 
and these electrons mostly distribute in the lower-energy 
regime compared to the high-energy regime overlapping 
the L valleys, the behavior of the SRT is mainly deter- 
mined by this part of electrons. In the approximation 
that all electrons are in the T valley, for E = 4 kV/cm, 
we have Tp sa 381 K according to the electron density at 
the peak, comparable to the corresponding hot-electron 
temperature T e m 659 K (T F /T e ~ 3/5). Meanwhile 
for E = 6 kV/cm, we have Tp w 563 K compared to 
T e fa 981 K at the peak (T F /T e ~ 3/5). We note that 
this effect of electric field on the density dependence of 
the SRT can be observed within current technology of 
optical orientation. 
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FIG. 4: (Color online) SRT r vs. electron density at tempera- 
ture T = 300 K for fixed electric field: E = 0, 4 and 6 kV/cm. 
Note that the scale of the case E = kV/cm is on the right 
hand side of the frame. 



C. Temperature dependence 

We now turn to investigate the temperature depen- 
dence of the SRT. In Fig. O the SRTs obtained are 
plotted as function of lattice temperature with electric 
field E = 0, 4 and 6 kV/cm and for three electron 
densities as in Sec. IIIA. In the field-free case, it has 
been formerly shown both experimentall y 10 i 18 ' 19 i 50 and 
theoreticall y 20 ' 51 that in n-type samples with low mo- 
bility, the SRT decreases monotonically with tempera- 
ture. From Fig. [SJa), it is seen that our result coincides 
with the previous ones. However, very different behav- 
iors are seen under high electric fields. Apart from the 
overall suppressed values of the SRT compared to the 
field-free condition, it is found that for electron density 
n e = 10 18 cm~ 3 , the decreasing rate becomes smaller for 
higher electric fields. Whereas for n e — 10 17 cm~ 3 , the 



SRT turns to increase with increasing temperature un- 
der electric field E = 6 kV/cm. The most interesting 
phenomenon is seen for the case with n e = 10 16 cm~ 3 . 
From Fig. Ela)-(c), one observes that the SRT decreases 
monotonically with increasing lattice temperature T in 
the absence of the electric field but increases monotoni- 
cally with it when electric field E = 6 kV/cm is applied. 
In between, for the case with E = 4 kV/cm, the SRT first 
increases and then decreases with increasing temperature 
T, with a peak at around T = 350 K. This peak in the 
temperature dependence of the SRT is very different from 
the one theoretically predicted by Zhou et al^ and ex- 
perimentally realized by Leyland et al.^ Ruan et al^ 
and Han et al£l for high-mobility samples in the field- 
free condition. There, the peak is solely caused by the 
electron-electron Coulomb scatterin g 38 ' 51 " — and appears 
in the crossover between the degenerate and nondegen- 
erate limits where the Fermi temperature Tp is compa- 
rable to the lattice temperature^ - — However, in n-type 
bulk materials, due to the strong electron-impurity scat- 
tering (hence low mobility), the Coulomb scattering is 
always less important and no peak is expected in the 
temperature dependence of the SRTJ£~— Moreover, the 
Fermi temperature corresponding to n e = 10 16 cm~ 3 is 
Tp s» 28 K, which is far below the lattice temperature, let 
alone the hot-electron temperature under the high elec- 
tric field. This further demonstrates the essential differ- 
ence of the peak observed here. 

This complex behavior of the SRT can be understood 
from the different lattice-temperature dependences of the 
hot-electron temperature at the different electron densi- 
ties and electric fields. Corresponding to Fig. [5^b) and 
(c), in Fig. O^d) and (e) we plot the hot-electron tem- 
peratures against the lattice temperature for E = 4 and 
6 kV/cm, respectively. Note that the same color and type 
of point are used for the corresponding electron density. 
By comparing the Tr-T curves in Fig. [5jb) with the r-T 
ones in Fig. \5^d) and those in Fig. [5fc) with Fig. 03e), 
direct correspondence of the temperature dependence of 
the SRT to that of the hot-electron temperature Tr is ob- 
served. The order of curves from top to bottom reverses 
in the corresponding two figures as higher hot-electron 
temperature indicates larger inhomogeneous broadening, 
and hence smaller SRT. Meanwhile, the decrease (in- 
crease) of the hot-electron temperature gives an increase 
(decrease) of the SRT in the lattice temperature depen- 
dence. Specifically, the monotonic decrease (increase) 
of the SRT with n e = 10 18 cm" 3 (n e = 10 16 cm" 3 ) 
in Fig. EJb) [Fig. 03c)] corresponds to the monotonic 
increase (decrease) of the hot-electron temperature in 
Fig. OJd) [Fig. Efe)]. Especially corresponding to the 
peak in Fig. G2h) for n e = 10 16 cm~ 3 , there is a valley 
in the lattice-temperature dependence of the hot-electron 
temperature in Fig. [5jd) . For the electric field- free case, 
the electron temperature equals the lattice temperature, 
so the behavior of the SRT in Fig. [5Ja) is also in the 
same trend as the high field cases. It is noted that in 
the discussion above, we focus on Tr instead of Tl. This 
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FIG. 5: (Color online) SRT r vs. lattice temperature T with 
(a) E = kV/cm; (b) E = 4 kV/cm and (c) E = 6 kV/cm. 
Hot-electron temperatures Tr and Ti, lis. lattice temperature 
T with (d) E = 4 kV/cm and (e) E = 6 kV/cm. The electron 
densities are n e = 10 , lO 1 ' and 10 10 cm in these figures. 



is because, as discussed in the previous subsection, the 
total spin relaxation is mainly determined by the F val- 
ley while most T- valley electrons stay in the lower-energy 
regime. Meanwhile, according to the the discussions in 
Sec. IIIA, in the variation of the hot-electron tempera- 
ture, the variation of inhomogeneous broadening is more 
profound than that of the momentum relaxation, thanks 
to the cubic k dependence of the T-valley Dresselhauss 
SOC i 43 ' 44 Therefore the variation of the inhomogeneous 
broadening mainly determines the variation of the SRT. 

We point out that the different behaviors of hot- 
electron temperature with the lattice temperature at 
different electron densities and electric fields are origi- 
nated from the temperature dependence of the energy- 
gain and loss rates of the electron system^ - — With the 
increase of lattice temperature, the mobility is reduced 
resulting from the enhancement of scattering, so does 
the the energy-gain rate; meanwhile the energy-loss rate 
decreases due to increase of phonon number^r— These 



two effects compete with each other and lead to the com- 
plex behavior of the hot-electron temperature in differ- 
ent conditions. A more detailed discussion is given in 
Appendix C. 



IV. SUMMARY 

In summary, we have investigated the multivalley spin 
relaxation in n-type bulk GaAs in the presence of high 
electric field by applying the fully microscopic KSBE ap- 
proach. The r and L valleys, which are relevant in de- 
termining the properties of spin dynamics for the high 
electric field applied in this work, are taken into account. 
The effect of the L valleys on spin relaxation is discussed 
and is shown to be very different from the quantum well 
system. We find that apart from the effect of directly 
manipulating the SRT, the high electric field can also 
effectively modulate the density and temperature depen- 
dences of the SRT. 

First, the SRT is found to decrease monotonically with 
the electric field for electron densities from the nondc- 
generate to degenerate limit. This monotonic field de- 
pendence of the SRT is very different from the previ- 
ous works in n-type GaAs quantum wells^— and is as- 
signed to the pronounced enhancement of inhomogeneous 
broadening from the field-induced drift and hot-electron 
effects thanks to the cubic form of the T-valley Dres- 
selhauss SOC in bulk. We show that, in despite of the 
very different strength of the SOC in different valleys, 
the evolutions of spin polarizations in the T and L val- 
leys share the same damping rate. This is demonstrated 
to come from the strong intervalley electron-phonon scat- 
tering and indicates the feasibility of exploring the prop- 
erties of the hot electrons in the L valleys through the 
phenomena detected in the T valley. Moreover, differing 
from the role of a "spin drain" of the total spin polar- 
ization in the two-dimensional system, we find that in 
bulk the L valleys serve as a "momentum damping area" 
where electrons are blocked from drifting to higher mo- 
mentum states due to the large effective mass. This tends 
to suppress the inhomogeneous broadening and in turn 
leads to a longer rather than shorter SRT compared to 
the case without the L valleys. 

As for the density dependence of the SRT, the for- 
merly predicted^ and experimentally observe d 13 ! 15 den- 
sity peak in the field-free condition is recovered and is 
found to be shifted to higher density regime with higher 
electric field. We attribute this to the electric-field- 
induced hot-electron effect. 

We also investigate the temperature dependence of the 
SRT in conditions with different electron densities and 
electric fields. The monotonic decrease of the SRT with 
increasing lattice temperature in the field-free condition 
coincides with the previous works. Nevertheless, the SRT 
is found to decrease more slowly with higher electric field, 
and even turn to increase monotonically with increasing 
lattice temperature in the condition with low electron 
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density and high electric field. More interestingly, a peak 
is predicted in the temperature dependence with low elec- 
tron density and moderate electric field, which is vastly 
different from the formerly discussed one in high mobility 
samples in the field-free condition. We point out that this 
peculiar behavior of the SRT originates from the temper- 
ature dependence of the energy-gain and loss rates of the 
electron system. 

Finally, we remark on the feasibility of the spin Gunn 
effect in n-type bulk GaAs. For the preferred electron 
density in Ref. 37 (n e = 10 18 cm -3 ), we note that the 
SRT is suppressed down to the value shorter than what 
required for the spontaneous spin amplification to appear 
under the electric fields where the charge Gunn effect 
appears. This fast damping of spin polarization overtakes 
the spontaneously generation process and makes the spin 
Gunn effect hardly to be realized in n-type bulk GaAs 
system. 



where 
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In these equations, p^ — pk and p^_ = 1 — pk- rii is 
the impurity density which equals the electron density 
in this paper and Zi = 1 is the charge number of the 
impurity. e[ r = fc£/(2m£) and e£* = k%J(2m* L ) + Etl 
with etl denoting the energy difference between the T 
and L points. N\y = [e (aAA '/ T) - 1] _1 is the Bose dis- 
tribution of phonons with frequency fl \\> . Nj\, = N\\> 
and N^ x , = 1 + N\\>. M^A'.q is the matrix element 
of the electron-phonon scattering with q standing for 
the phonon wave vector. Here we take into account 
the intra- and intervalley electron-longitudinal-optical 
(LO) phonon scattering in and between the T and L val- 
leys, respectively! 35 ' 41 ' 59 " — For the intravalley electron- 
phonon scattering, we have M 2 T r = e ilrr [ fl ' x 1> K ° - 1 and 
n 2 



2eoq 2 



Also we have M 2 L . 



M ir, q 



2d£2 r 



for the T-L intervalley electron-phonon scatter- 

ing and M^, L . q = 2dn ' L jf for L-L intervalley electron- 
phonon scattering. V q is the screened Coulomb potential 
under the random phase approximation^ 3 - 



Appendix A: Scattering Terms in KSBEs 



The scattering term d t p\k x l sca t includes the contribu- 
tions from the electron-impurity scattering dtpxk x | e i! the 
electron-phonon scattering d t p\k x \ cp and the electron- 
electron Coulomb scattering d t p\k x \ ee 



dtPXk x | scat = d t p\k x | el + d tP \k x | ep + d t p\k x | ec , (Al) 



V a 



(0) 



l-V/ q (0) PW(q) 



where 



P«(q) = 



E /A(k A +q).cr fxk x ,a 
£ X A : 

A,k A: er k A +q k A 



(A5) 



(A6) 



with Vq * 1 — e 2 /( e o K o<? 2 ) denoting the bare Coulomb po- 
tential and /Ak,o- being the electron distribution of spin-cr 
band. It is noted that in Eq. (|A4I) . we include both the 
intra- and intervalley electron-electron Coulomb scatter- 
ing. All parameters appear in these equations are listed 
in Table. I. 
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FIG. 6: (Color online) Electric field dependences of (a) drift 
velocity v\ for electron density n e = 10 16 cm -3 , (b) mobility 
[i, (c) ratio of electron densities in the L and F valleys and 
(d) hot-electron temperature for three electron densities n e — 
10 16 , 10 17 and 10 18 cm" 3 . The lattice temperature is T = 
300 K and the experimental data are taken from Ref. 1.321 . 



Appendix B: Drift velocity, Mobility, L- valley 
Occupation and Hot-electron Temperature under 
Electric Field 



The hot-electron temperature T e is obtained by fitting 
the calculated steady-state electron distribution of each 
valley with the drifted Fermi distribution /(ek A ,B) = 
[exp((e kAvE - n\)/T\) + l] -1 . Here e kA ,£; = (k A - 
k°) 2 /2m A is the energy spectrum shifted by the elec- 
tric field with being the drift momentum and [i\ is 
a fitting parameter denoting the chemical potential in A 
valley. It is found that electrons in the T valley carry 
two temperatures, one in the higher-energy regime that 
overlaps with the L valleys (labeled as Tj,) and the other 
in the lower-energy regime (labeled as Tr). Electrons in 
the L valleys share the same temperature as those in the 
higher-energy regime of the T valley due to the rapid ex- 
change of electrons thanks to the strong intervalley scat- 
tering. In Fig. |5[d), we plot Tr and Tl against the elec- 
tric field at lattice temperature T = 300 K. It is seen that 
electrons in the T valley are easier to be heated due to 
the smaller effective mass. Moreover, by increasing the 
electron density from 10 16 to 10 18 cm -3 , Tr is effectively 
reduced while Tj, stays almost unchanged. The underly- 
ing physics is that by increasing the electron density, the 
electron-impurity scattering is enhanced. This tends to 
reduce the drift velocity in each valley and suppress the 
ability of electrons to gain energy from the electric field, 
and therefore reduce the hot-electron temperature^^— 
Since the electron-impurity scattering is the leading scat- 
tering in the T valley^ whereas the intervalley electron- 
phonon scattering is dominant in the L valleys^ 8 - the drift 
velocity, and hence also the hot-electron temperature of 
the L valleys, are less affected by the electron density 
compared to those of the T valley. 

Appendix C: Energy-gain and Loss Rates 



In order to have an overview of the electric proper- 
ties influenced by the high electric field, the steady-state 
drift velocity V\ of each valley as well as the total drift 
velocity are calculated by varying the electric field from 
to 8 kV/cm. In Fig. UJa), we plot the drift velocity 
v\ as function of electric field with n e = 10 16 cm -3 and 
T = 300 K. The negative differential electric conductance 
can be seen from the field dependence of the total drift 
velocity and good agreement is reached with the experi- 
mental results. 32 We also calculate the field dependences 
of the mobility and the high-valley electron population 
for three different doping densities. From Fig. EJb), one 
finds that the mobilities first decrease slowly with the in- 
crease of electric field and then more rapidly when the 
electric field is increased over E — 4 kV/cm. This can be 
understood with the help of Fig. |6jc) where the ratio of 
electron densities in the T and L valleys is plotted against 
the electric field E. It is seen that the population of elec- 
trons in the L valleys is negligible when E < 2 kV/cm 
and approaches 10% when E ~ 4 kV/cm. This con- 
tributes to the faster decrease of mobilities in Fig. HJb) 
and leads to the negative differential electric conductance 
in Fig. HJa). 



For a semiconductor system under uniform electric 
field, the electrons accelerate before they are scattered 
and thus gain energy from the electric field. Mean- 
while due to the electron-phonon scattering, the elec- 
trons transfer energy to the phonon system. In steady 
state, the electron energy-gain rate equals the energy- 
loss one *^r— 

We calculate the energy-gain and loss rates in n-type 
bulk GaAs by including only the T valley. The energy- 
gain rate rj (in unit volume here and hereafter) reads^ 

77 = en e fiE 2 (CI) 

and the energy-loss rate w is given byiil 

u = 2^r> q , x |M(q, X )| 2 n 2 (q,r! q , x + W o) 



q,\ 



with 



ra (gsbX)_ n ( "q.x + "o , 
T v T„ 



n 2 (q ;W ) = 27r£[/(e k) T e )-/(e k+qj T )] 

k 

X <5(£k+q - £ k + W). 



(C2) 



(C3) 
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In these equations, ^ q , x is the phonon energy with mo- 
mentum q and mode \- Note that here we only need 
to take account of the intravalley electron-LO phonon 
scattering with |Mrr,q| 2 given the Appendix A. ljq = 
q • Vd with Vd denoting the drift velocity. n{x) = 



1 / (e x — 1) stands for the Bose distribution and f(x, T e ) = 
l/[e^ x ~^' T ' + 1] is the Fermi distribution with T e and fj, 
being the electron temperature and the chemical poten- 
tial, respectively. 
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FIG. 7: (Color online) The difference of energy-gain and loss rates u — r\ against the lattice temperature and hot-electron 
temperature, (a) n e = I0 16 cm -3 and E = 4 kV/cm; (b) n e = 10 18 cm" 3 and E = 4 kV/cm; (c) n e = I0 16 cm -3 and 
E — 5 kV/cm. The dashed curves indicate the points where u) = r], hence are just the plots of T e -T in equilibrium. 



From Eqs. (|C1|) and (|C2|) . we calculate the energy-gain 
and loss rates in three cases: (i) n e — 10 16 cm -3 with 
E = 4 kV/cm; (ii) n e = 10 18 cm" 3 with E = 4 kV/cm 
and (iii) n e = 10 16 cm" 3 with E = 5 kV/cm. The lat- 
tice temperature is varied in the range T = 200 ~ 500 K 
and the hot-electron temperature T e in the proper range 
according to the electron density and the electric field. 
The drift velocity used in the calculation are obtained by 
solving the KSBEs. In Fig. [7J we plot ui — r\ against T 
and T e for the three cases. The points where tu — r; = 
are indicated with dashed curve, which is exactly the hot- 
electron temperature T e versus to the lattice temperature 
T in the steady state. By comparing the dashed curves 
in Fig. [TJa) and (b) with the corresponding curves in 
Fig- Eld), one notices that qualitatively, good agreement 
is reached. We note that for case (iii), due to the large 
electric field, this model does not hold so well in describ- 
ing the genuine system since the L valleys start to play 
an important role and serve as the "momentum damp- 
ing area" where electrons are hardly drifted. Therefore in 
the simplified model, the drift effect is overestimated and 
the intervalley electron-phonon scatterings, which serve 
as additional energy-loss channels, are missing. These 
lead to the overestimation of T e . However, the qualita- 
tive behavior of hot-electron temperature with the lattice 
temperature in the presence of extremely high electric 
field is still captured: T e decreases with increasing T. 

The different behaviors of T e with T can be under- 
stood as follows. The increase of the lattice tempera- 



ture induces two main effects: (1) It reduces the mobil- 
ity by enhancing the electron-phonon scattering, which 
in turn leads to the decrease of the energy-gain rate [see 
Eq. (|C1[) ]. This tends to reduce the hot-electron temper- 
ature. (2) It also modulates the relative importance of 
phonon-emitting and absorbing processes in the electron- 
phonon scattering and reduces the energy-loss rate for the 
electron system (by reducing the temperature difference 
between the electron and phonon systems). This tends 
to "heat" the electron system. The competing of these 
two factors contributes to the complex T-dependence of 
T e , hence also that of the SRT. In Fig. [71 we have shown 
three typical cases investigated. In Fig. [TJc) with high 
electric field and low electron density, the hot-electron 
temperature T e is high above the lattice temperature T. 
As a result, the small increase in T does not affect the 
energy-loss rate much but does enhance the scattering. 
This results in the leading role of the reduction of the 
energy-gain rate, which in turn leads to the decrease of 
T e . Whereas in Fig. UJb) where T e is much closer to 
T compared to that in Fig. [TJc), by increasing the lat- 
tice temperature, the temperature difference between the 
electron and phonon systems is effectively reduced and 
hence the "heating" effect is more efficient. This leads 
to the increase of the electron temperature. In between, 
for the case with relatively low electric field and low elec- 
tron density, a nonmonotonic temperature dependence is 
expected, which is just the case in Fig. UJa). 
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